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The recent introduction of a deformed non-minimal version of the noncommutative Standard 
Model in the enveloping-algebra approach, having a one-loop renormalisable gauge sector involving a 
higher order gauge term, motivates us to consider the possibility of extending the fermion sector with 
additional deformations, i.e. higher order fermionic terms. Since the renormalisability properties of 
the fermion sector of the model are not yet fully known, we work with an effective fermion lagrangian 
which includes noncommutative higher order terms involving a contraction with the noncommutative 
8 tensor aside from the star products, so that these terms annihilate in the commutative limit. Some 
of these terms violate CPT in the weak sector, and some violate CP in the strong and hypercharge 
sectors. We apply this framework to the reevaluation of the decay rates of quarkonia (qqi = J/tp, T) 
into two photons. These decays, which are forbidden in the ordinary Standard Model, had been 
previously studied as possible signals for noncommutativity, but not in the framework of the better 
behaved deformed non-minimal version of the noncommutative Standard Model. Weak CPT or 
strong-hypercharge CP violating interactions do not contribute to the result. If the parameters of 
the model take natural values, for the vast majority of configurations the resulting branching ratios 
are enhanced with respect to their values in the minimal version of the noncommutative Standard 
Model. Also for more than half of the parameter space, the rates are larger than the maximal 
rates that were calculated in the undeformed version of the non-minimal noncommutative Standard 
Model. Tuning the dimensionless parameters the predicted branching ratios can fall within the 
current experimental bounds. 

PACS numbers: ll.10.Nx, 12.38.-t, 12.39Dc, 12.39.-x, 14.20-c 

INTRODUCTION 

The Standard Model (SM) of particle physics and the theory of gravity describe very well, as far as we know today, 
all physical phenomena from cosmological processes to the properties of subnuclear structures. Nevertheless, at 
extreme energies and/or very short distances -at the Planck scale- this theories fail to be compatible, which motivates 
the study of modified or alternative space-time structures that could help to solve the above mentioned difficulties 
or at least shed some light on them. These modified space-time structures arise in such settings as the quantised 
coordinates in string theory or in the general framework of deformation quantisation. The idea of noncommutative 
(NC) space-time, which can be realised in both of the above settings, has recently found more and more interest. In 
this paper we deal with noncommutative theories defined by means of the enveloping algebra approach, which allows 
to define gauge theories with arbitrary gauge groups, in particular that of the Standard Model. The research on these 
theories so far has successfully dealt with some theoretical and also phenomenological aspects, which might allow the 
confrontation of the theory with experiments. 

On the theoretical side, in the enveloping algebra approach developed in [1] following the ideas of the seminal paper 
[2] , we can emphasise the construction of noncommutative minimal and non- minimal versions of the Standard Model 
(mNCSM and nmNCSM) [3-5], and GUTS [6]. Of particular interest is the deformed version of the nmNCSM with 
renormalisable pure gauge interactions discussed in refs. [7, 8], which involves an extension of the pure noncommutative 
gauge lagrangian F * F with a deformation or higher order term. These noncommutative extensions of the SM are 
anomaly free [9, 10]. In contrast with the gauge sector, the renormalisability of the fermion sector has not yet 
been completely addressed, despite some encouraging partial results [11]. Nevertheless, the results in the pure gauge 
sector motivate to consider the possible effects of introducing higher order fermionic terms; in this paper we will 
initiate the study of an extended fermion sector by introducing a general class of the fermionic deformations in the 
noncommutative action of the nmNCSM and obtaining some of the corresponding Feynman rules. 

On the other hand, on the phenomenological side, it is known that noncommutative field theories can predict non- 
zero rates for processes which are forbidden in the Standard Model due to the Lorentz invariance and Bose symmetry 
(Landau-Pomcranchuk-Yang or LPY theorem). These new effects follow from the violation of Lorentz invariance in 
the presence of noncommutativity, and the observation of these forbidden decays could be taken as a signal for it. 
Some of the SM forbidden processes that have been studied, at the level of tree diagrams, include Z — > "f"f/gg [4]; 
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J/\I/,T — > 77 [12] and if — > 7T7 [13]. For other phenomenological studies, including limits on the noncommutative 
scale AnCj see [14]. The recent introduction of the deformed nmNCSM calls for a reevaluation of the results, which so 
far has only been done for the Z — > 77 process [15]. In this paper we will apply our framework of a deformed fermionic 
sector to study the quarkonia decays into two photons J r /\P ', T — > 77. There are other processes involving quarkonia 
in which noncommutativity induces rare decays. An example is the decay into two photons of quarkonia polarised 
in some direction, say the x 3 axis. In such polarised rate the contributions proportional to the third components 
(Eg) 2 and (Bg) 2 should be enhanced by a large factor, similarly as in the Z — > 77 case [15]. In addition, one could 
study decays of quarkonia in other gauge bosons. In general, we expect that any quark-antiquark state with the 
same quantum numbers as those of J/^f and T (I G (J PC ) = 0~(1 )) will decay into two gauge bosons through 
noncommutative interactions. We restrict ourselves to quarkonia decays into two photons because their detection 
could in principle be achieved within very high-resolution calorimeters by applying more stringent conditions in the 
selection of the photon candidates when searching for 77 events. We could also have decays into a pair of gluons, but 
these would hadronise into hadron jets, whose detection would be much more problematic due to a lack of localisation 
of the jets or interference with other signals. Of course there could be also decays into two Z bosons, but, on the one 
hand, they would not be produced for quarkonia at rest due to the heaviness of the Z bosons, and, on the other, these 
would rapidly decay into other particles. 

The paper is organised as follows. First, we give some more detailed fundamental and phenomenological motivations 
for our work, and we introduce the theoretical framework. Next we introduce the action of the renormalisable pure 
gauge sector of deformed nmNCSM, after which we deal with the matter sector, extending the fermion action with 
higher order deformations and examining their C,P,T transformation properties. Following this, we study whether 
field redefinitions in the fermion sector can help to get rid of some of the extra terms introduced. Finally we give the 
Feynman rules relevant to the calculation of the quarkonia decay rates and we present our results, which are then 
discussed. 

MOTIVATION AND FRAMEWORK 

Our main goal in this paper is, following the recent introduction of the deformed version of the gauge sector of 
the noncommutative Standard Model (NCSM), to consistently define the action of an extended fermion sector of the 
nmNCSM from an effective-theory point of view, by adding higher order terms, and derive relevant Feynman rules. 
This is certainly important for future investigations about the renormalisability properties of entire NCSM, which are 
interesting in their own right and could introduce constraints on our "effective" deformed fermion sector. The second 
goal of this paper is to apply this framework to the calculation of the decay rates of the simplest processes which are 
forbidden in the ordinary SM, that is the C symmetry violating decay of quarkonia into two photons, (qqi = J/ip, T). 
These processes are of theoretical interest because their tree-level contributions come from two different types of 
diagrams, these being s-channel gauge bosons exchanges -involving interactions from the pure gauge sector of the 
NCSM- and t-channel quark exchanges -involving interactions from the fermion sector. Of course, it should not be 
forgotten that this proposed processes are also important because they violate the LPY theorem that holds in the 
ordinary SM, and thus their hypothetical detection could be taken as a possible signal of noncommutativity. 

The first construction of the NCSM was undertaken in ref. [3], where it was already noted that there was an 
ambiguity in the choice of traces for the gauge kinetic terms, leading to a minimal version (with traces taken in the 
adjoint representation), and a non-minimal version including traces over the representations of all the massive particle 
multiplets charged under any of the gauge groups [4] . The interaction vertices and Feynman rules of these models 
were further analysed in refs. [5]. As was said, some of the processes studied in this framework include Z — > 77/gg 
[4, 15]; J/*, T 77 [12], K -> Try [13] and [16]. 

At one loop order, the investigations so far were concerned with renormalisability properties [7, 8, 15]; remarkably, 
it was found that the pure gauge interactions of the nmNCSM could be rendered one-loop renormalisable at first order 
in the noncommutativity parameters 6^ u by adding an extra deformation term to the lagrangian involving only gauge 
fields and their derivatives contracted with one 6^ v . This term was introduced at the level of the noncommutative 
action [8, 15], yielding an extended version of the nmNCSM [7]. This result was reached by considering only gauge 
field contributions to the loop integrals. When matter fields are included, the results of ref. [17] show that, when 
computing the matter contributions to the one-loop diagrams with external gauge fields in a generic noncommutative 
gauge theory with Dirac fermions or complex scalars, the divergences can be absorbed in the bare lagrangian whenever 
the representations of the matter fields are included in the choice of representations for the traces in the pure gauge 
terms of the action. 

This result motivates the following comments. First, it suggests that there is virtually no hope of getting renormal- 
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isability of the gauge sector in the mNCSM case when Dirac fermions run in the loops, since the model only involves 
traces in the adjoint representation of the gauge fields. However, in the pure gauge sector of the nmNCSM [4, 5] the 
traces include the representations of the matter fields, so that one could say that it is likely that the nmNCSM gauge 
sector will still be renormalisable if one includes the effects of the matter fields in the loops. This motivates even more 
to favour the use of the nmNCSM, and in particular its deformed version, over the mNCSM. 

Regarding the renormalisability properties of the nmNCSM, we would like to make some clarifying observation. 
The NC SU(N) pure gauge theory [8] showed renormalisability for two choices of the free deformation parameter a, 
associated with the higher order term, S^(a), that is present in the deformed nmNCSM: a = 1 and a = 3. However, 
this result cannot be extrapolated directly to the NCSM, because the gauge fields mix after the SW map and it is not 
a sum of NC SU(N) theories. In fact, the gauge sector of the nmNCSM is only renormalisable for a — 3, even more, 
the one-loop quantum corrections are finite for a = 3. 

The above arguments clearly favours the nmNCSM over the mNCSM if one wants to have a better behaved theory, 
where the quantum corrections are more under control even when working in an effective theory approach. In 
particular, predictions of the theory involving pure gauge boson interactions will be more robust under changes of 
scale. 

All these new results call for a reevaluation, in the framework of the deformed nmNCSM, of the possible signals 
for noncommutativity that were commented upon before. As was said before, this has already been done in ref. [15] 
for the Z — > 77 decay. In this paper we will focus on the disintegration J/\P,T — ► 77. Now, since the renormalis- 
ability properties of the fcrmion sector of deformed nmNCSM are essentially unknown, despite promising results for 
some diagrams involving chiral fermions [11], and since the results for the pure gauge sector make us expect that a 
renormalisable fermion sector could only be achieved by adding deformed fermionic contributions to the action, it 
makes sense to treat this sector effectively and consider, aside from the fermion lagrangian employed in both mNCSM 
and nmNCSM, higher order contributions involving contractions with 9^ outside the star product and compatible 
with the noncommutative gauge symmetry. We will demand that these deformed fermionic contributions do not 
introduce further violations of parity than those coming from the ordinary SU(2)l gauge fields; also, we will require 
that they do not alter the tree level fermion 2-point function. This modification of the propagator could be achieved, 
for example, with terms like iTp6 a P^p{D 2 , D a } *ip. The problem is that, when computing S-matrix elements, one 
ordinarily uses the LSZ formalism, which is constructed by using Lorentz invariance and implies that the S-matrix 
elements are obtained by identifying the poles in momentum space of some Green functions. This happen at points 
satisfying the Lorentz invariant constraint p 2 = m 2 for some m. In noncommutative spacetime Lorentz invariance is 
broken, which means that the poles of the Green functions need not satisfy the constraint p 2 = m 2 ; a ^-dependence 
may appear and then one should be more careful when computing S-matrix elements. An alternative definition is 
needed, perhaps along the lines of the work done in ref. [18] for noncommutative theories formulated without SW 
maps. Nevertheless, if the tree-level propagator still has a pole of the ordinary type, we expect that the usual way 
of deriving matrix elements will be valid to some approximation. Here we pretend to compute S-matrix elements at 
tree level only. 



PURE GAUGE SECTOR 
Renormalisable pure gauge sector action 

The pure gauge part of the deformed nmNCSM is given by [15]: 
S g = Sf n + S»{a), 

Sf n = ~\^ j d 4 xF„ v * Sf(a) = ^Tr J d l xhff iV * % v * F pa * (1) 

where the trace Tr is taken over all the particle representations. As usual, the Moyal-Weyl ^-product is given by 
f * g = f (exp I hO^ d ^ d g, which implements the space-time noncommutativity as x^ * x v — x v * x^ = ih6 tiv . 

The noncommutative deformation parameter h — 1/A^ C sets the noncommutative scale. The noncommutative field 
strength F^ depends on the enveloping- algebra valued noncommutative gauge field as 



% v {x) = dX - d v % - i[% * V v ], 



(2) 
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and in turn, V p depends on the ordinary gauge bosons through the Seiberg-Witten map 

%{x) = V^x)-j6°r{V a (x),d p V^x)+F ptl (x)} + 0(h 2 ), (3) 

3 8 

V„{x) = g'A p (x)Y + g J2B f >, a (x) T L+9sJ2 G ^W T s> ( 4 ) 

a=l 6=1 

where V^(x) is the Standard Model gauge potential taking values in the Lie algebra of SU(3)c x SU(2)l x U(1)y- 
The pure gauge action defined by eq. (1), after expanding to order h with the SW map of eq. (3) leads to 



S g = Tr J d 4 x (-^F^F^ + h6» u (| F pv F pa - F Mp F vtT )F pc 



(5) 



After taking traces in the above action over all massive particle representations with different quantum numbers that 
appear in the total lagrangian of the model with covariant derivatives acting on them, the gauge action (5) produces 
triple neutral gauge boson interactions [4, 5] which are not present in the mNCSM; in this paper we are interested in 
the 777 and Z77 couplings, which arise from the following terms in the lagrangian [15] 

e 



-nmNCSM / 



sin2<V K^hO^A^ {aA^Apv - AA PP A UIJ ) 



C 



^ CSM (a) = -sm26 w K z ^h6P° 2Z^(2A pp A va - aA pv A pa ) + 8Z^A^ A va - aZ pa A^A^ 



(6) 



See fig. 1 for details on the allowed values for the constants i^ 777 and if z 77 and their dependence on other parameters 
of the model [4]. The remarkable result of refs. [7] is that the gauge action (5) is one-loop renormalisable up to order 
h for a = 3, with the 0{h) quantum corrections being finite. Thus we will be mainly interested in the a = 3 case. 




-0.5 -0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 

FIG. 1: Allowed values for the couplings iCz 77 and if 777 in the nmNCSM. 

MATTER SECTOR 
Minimal fermion action 



Next we turn our attention to the fermion sector. In both the mNCSM and nmNCSM models, the chosen minimal 
fermion action was 



d Xtp * ilp * 1p + SVukawa, 0^ = 0^- iV p *, 



(7) 
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which depends on the ordinary fields V p and ip through the SW map of eq. (3) and the fermion SW map 

i> = 1> - \ ap (v a d p - % - [V a , Vp]^ 4, + 0{h 2 ) . (8) 

The Yukawa sector in the NCSM, when expanded in terms of ordinary fields, involves quite complicated interactions 
[5]. However, if we consider the noncommutative contributions to the fermion field interactions in the QED sector 
only, (as turns out to be sufficient for our purpose of obtaining rate for the quarkonia decay into two photons), we 
can use, as noted in [5] and used in [12], the simplified fermion Lagrangian 

S^a ~ fat* (ip A - m,) *£ = Sffi + S^ A + 0(h% 

S^l = J d^(i]/) A - m s ) V, Sl, A = ~J d 4 x^A^ (&»"> D A - m f 6^, (9) 
where D A = d p — ieA^ and 



Deformed fermion action 



Motivated by the unknown rcnormalisability properties of the NC fermion sector, and since in the case of the 
gauge sector the addition of the higher order term dependent on the free deformation parameter a made the model 
surprisingly well-behaved, we will consider adding deformation terms to the minimal fermion lagrangian (9), either to 
treat it in an effective theory approach or to prepare the grounds for future investigations of renormalisability. Thus, 
we set to find all the possible deformations contributing to the fermion action, and satisfying the following conditions: 

(a) They are real and include two fermions fields, 

(b) are invariant under noncommutative gauge transformations and thus involve star products, noncommutative gauge 
covariant derivatives and noncommutative field strengths, 

(c) involve a contraction with a tensor outside the star product — as in the a-dependent term in cq. (1) — , 

(d) they do not alter the tree-level 2-point function of fermion propagator, 

(e) they only generate P-symmetry violating contributions coming from the ordinary weak SU(2) gauge fields, 

(f) include zero or positive powers of some mass parameter. 

In order to write down all the possible terms satisfying the conditions above, we have to consider a basis of matrices 
in spinor space. Since the 7 matrices are constrained by the identity {7^, 7^} = 2g fJiU , we can consider a basis formed 
by antisymmetrised products of 7 matrices: {/, 7^, <7 M!y , 7 A " yp , 7 5 }, where we are using the following definitions: 

o»» = I(yy_yy), 

yii/p _ I^yxyy> _|_ yyy _|_ yyy _ yyy _ yyy _ yyy^ 
yyy = ^y _ g ^ + «^y _ le ^^ la , 

i 5 = -JwtYiY. (11) 

With this in mind, it can be seen that the terms satisfying the conditions stated above are given by sums of integrals 
of the following monomials U multiplied by real coefficients: 



tl 




ti 






*3 


= th6^hl p (2F aP *D P + {V p F a p))*4,, 


U 


= ihe a ^ * Y(2F l3p *D a + (t> a F Pp )) * 




t 5 


= ^1*7^*^ + (T^F^))*^, 


t e 






*7 


= ihd a ^* 1Si ^{V a F pv )*i>, 


t& 


= he ati ^^ lap f(2F lip * i> + {iPFw)) 


■kip, 


*9 


= h6 a ^* lp ^{2F plJ *b a + { r b a F pv )) -k "0, 


tio 


= he a ^* lf3 ^(2F va *D^ + {T) p F va ) ) 


■kip, 


til 


= mh6 a0 $*F al3 *$, 


tl2 


= imh9 a ^*a[3 p F ap *$, 




tl3 
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with 9 a/3 = ±e a/3pc '6 ' pcr , and = — iV^-k , 2? M = — z[V^ * ] . Note that we can identify the mass parameter m 
with any of the fermion masses to/; this means no loss of generalisation because in principle we will be adding these 
mass-dependent terms to the action multiplied by arbitrary dimensionless coefficients. 

In order to check the C,P,T transformations of the above terms, we consider the following: 

• d^ v transforms under discrete space-time symmetries as a U(l) field strength F^ v [6]. 

• The ordinary fields are the ones that define the theory, so that we have to deal with their C,P,T transformation 
properties. However, the analysis can be simplified because of the fact that the SW maps that we use are such 
that the C,P,T transformations of the noncommutative fields are equal to their commutative counterparts [6]. 
However, there is a subtlety since the noncommutative vector field in the matter representations, given the 
expansion of eq. (4), includes some chiral projectors that come together with the weak gauge fields B Mjtl (x). 

Taking this into account, it can be seen that the above terms are P invariant save for the contributions involving the 
ordinary SU(2) gauge fields, as was required from the start. Moreover, all terms in eqs. (12) are CPT invariant except 
for in,ii2,£i3- In the terms in, £12,^13 the CPT violating contributions come exclusively from the chiral projectors 
of the weak gauge fields. However, in the strong and hypercharge sectors the contributions from all terms in (12) 
remain CPT invariant. Interestingly, the terms ti, ti, t$, tg, tio, £12 originate both C and T violations in the strong 
and hypercharge sectors, with CT conserved in these sectors. All other violations of C, P or T come from the weak 
fields exclusively. We will not worry about these violations of discrete symmetries and will proceed considering all of 
the terms in eq. (12). Thus, we will consider an additional piece of the action given by a sum of the above terms, 

S"(xi)= [d 4 xJ2xiU, x,eR, (13) 

i 

where Xi's are the fermion sector free deformation parameters, to be constrained via considerations of renormalisability 
or by phenomenology. 

Before plunging into the computation of the Feynman rules and their phenomcnological application, we would like 
to analyse whether any of the terms in eq. (12) can be reabsorbed by means of field redefinitions. This is interesting 
for two reasons: 

- Possible future investigations about renormalisability would require to take into account the effects of field redefini- 
tions in the action, which we will compute here in the fermion sector. 

- Field redefinitions should not affect S-matrix elements, and since we aim to calculate quarkonia decay rates it is 
desirable to eliminate as many of the Xi parameters in eq. (13) as possible. 

To proceed, we consider the most general field redefinitions of order 9 of the commutative gauge and fermion fields, 
involving zero or positive powers of the mass parameter and not affecting the tree-level propagator, and compute the 
associated change in the action. We parameterise the most general field redefinitions up to order as 

5v tl = ^y i 5 i v ll , 5i) = ^2z i 5 i i), (14) 
where yi are real, Zi complex, and 5 t v li and 5 l ip are given next 

6\ = e a ^F a0 , £ 2 v lt = e li a v v F va , 

,5V = e a ^F a ^, ^ = e^a "F ap i,, 6 3 ^ = e a ^ 5 F a ^. (15) 

Since the field redefinitions are of order h, to see their effect on the action to this same order we only need to compute 
the variation of the O(h ) action, which for fermions and vector fields is just 



qSM 



S s g M + J d 4 xiP (ip - to/) i> + S Yukawa ; (16) 

we will ignore the Yukawa couplings since they do not contribute to the tree-level quarkonia decay amplitudes that 
we want to compute. The change of the total action turns out to be 

5 (s% M + J d 4 x^{ifl - to/) ^ = Jd^Y^iVt^C + Rez, 5 Zi ' R C + Imz, 6 Z > J C) + 0(h 2 ), (17) 
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with 



l C = t u 8 V2 C = -t 2 , S Zl ' R £ = t 3 -2t 



S' 1 ' I C = -t 1 , S Z2 ' R £ = -t 4 + t B - ^7, <5 Z2j £ = -^i-t 2 -ii -2i 12 , 

^3,i?£ = _ ^3,/£ = _^ _ tj _ 2tlZi (18) 

where the tj's are given in eq. (12). In the language of eq. (13), the above result is equivalent to the following change 
in the action (xi): 



S"{xi)^S^{xi+5xi), (19) 



with 



Sx\ 


= 2/1 - Imzi 




<5x 2 


= -y 2 - Imz 2 , 


8x3 


= Rezi , 


8x4 


= -Rcz 2 , 




8x5 


= Rcz 2 , 


8X6 


= -lmz3 , 


8x7 


= -^Rez 2 - 


Imz 3 , 


8x s 


= -Rez 3 , 


8xg 


= -Rez 3 , 


8x w 


= -Imz 2 , 




Sxu 


= -2Rezi, 


8X12 


= — 2Imz 2 


8x 13 


= — 2Imz 3 . 













(20) 

In order to check whether the field redefinitions above can be used to eliminate some of the terms U of the action 
in eq. (13), we have to examine the system of equations that follows 

8xi[yj,Rezj,Imzj] = ~x l , (21) 

where Xj are to be treated as fixed and z/i, Rez^, Imzi as the unknown variables. Now, the 13 x 8 matrix associated 
to the previous linear system of equations can be seen to be of rank 7. This means that at most 7 of the U terms in 
the action can be eliminated with adequate field redefinitions, with 6 terms remaining. It is easily seen that a viable 
choice of 6 terms that survive the field redefinitions is given by t4,t$,t§,ts,tiQ,t\i note that, as was commented 
before, tu violates CPT in the weak sector and is>*io violate C and T in the strong and hypercharge sectors. 

FEYNMAN RULES TO COMPUTE S-MATRIX ELEMENTS 

From the discussions above it follows that, in order to compute the amplitudes of the desired physical processes, 
the relevant pieces of the action are those given by the ordinary Standard Model fcrmion-photon and fermion-Z boson 
interactions, plus the following noncommutative ^-dependent interactions: the three boson interactions of eq. (6), the 
fermion-photon interaction in eq. (9), and the contribution to the action of the terms of eq. (13) surviving the fermion 
field redefinitions, i.e., 



$(xi) = Jd 4 x x 4 t 4 + x b t b + x%t% + x 8 t 8 + xioho + £11*11 . (22) 



The noncommutative interactions originate, among other vertices, a triple neutral gauge boson vertex, a 2 fermion- 
photon vertex and a 2 fermion-2 boson vertex. From the modified gauge and fermion actions, (1,5,6,9,22), 

S = S g + = S™ m + S»(a) + Sf n + S$ ( Xi ) + 0(h 2 ), (23) 

we obtain the following Feynman rules: 




■iee 



7 M - i\ K(6% P P P + 0%" k P - (1 + 4z n ) m f 6*" ) 



+ 8ix 8 g^6^^ v + %ix w 6™ la «*, 
~ d H P = \(°h P - ^ VP ) - W*"^ e q = |, -1 for quarks, (24) 
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A ^ fci k 2 K 




~e 2 el [6?r] af) (h-k 2 y. 



fcgH A p 



2esm(29w)K 711 h Q% up [a; k\,k 2 , k 3 



(25) 



(26) 



fcsH A p 



-2e sm{26 w )K Zll h Q 3 lvp [a; k u k 2 , k 3 ] . 



(27) 



Here, 9 pup was given in eq. (10), while the tensor <d 3 iup [a; k\, k 2 , k 3 ], is given by [15]: 

0f [a; k u k 2 , k 3 ] = - (h6k 2 ) [(h - feOV" + (fc 2 - hTg»P + (fc 3 - fci) VI 

- 0"" [k p (k 2 k 3 ) - k p (fcifcs)] - 0"' [A# (fc 3 fci) - ^3 ( fc 2fci)] - W (fcife) - (*3fe)] 

"^1 



(28) 



+ [g vp kj - k»k p ] + (9k 3 y [g vp k\ - k v 2 k 

+ (9k 3 y [g» p k\ - k p k{\ + (6hy [g pp k 2 3 - k»k p ] 

+ (ekt) p [g pv k\ - k p k v 2 \ + {6k 2 ) p [gT k\ - k p k\\ 

+ (a- l)((0fci)" [ff"" (fc 3 fc 2 ) - fc£*£] + (6k 2 y [g^ihh) - k^k p ] + (6k 3 ) p [^"(feto) - k p k v ] ) 



APPLICATION OF THE PROPOSED FRAMEWORK TO QUARKONIA DECAYS: q qi -» 77 

Quarkonia decay amplitudes 

The diagrams that contribute to the quarkonia decay amplitude are shown in figures 2 and 3. Note that the vertices 
with fermions including noncommutative effects (with black dots), only involve photons; that's why we may safely 
use (9). The diagram in fig. 3 also involves a fermion-fermion-Z boson vertex, but it is given by the ordinary SM 
contribution. 

Using the Feynman rules (24-27) obtained from the modified gauge (6) and fermion (23) actions, we have evaluated 
the diagrams from Figs. 2 and 3, yielding amplitudes that we call Ai and A 2 , respectively. 

As in ref. [12], in order to hadronise the free quarks into the quarkonium bound state, we apply the following 
prescription for the transition amplitude of the operator qflf- (q = c,b and i,j are colour indices) from the vacuum 
to the quarkonium state: 



(o\Q?^\m(P)) = --^=- W + mW p <% , (29) 

where [^^(O)! represents the quarkonia wave function at the origin defined in [12], 
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We use a collinear approximation for the quarks in the quarkonium state, and we identify the mass of the quarkonium 



q e 



e(P) 



e{P) 



e(fci) 

e(fc 2 ) 
e(*i) 
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e(P) 




FIG. 2: Contributions to the Ai(qqi —> 77) amplitude. 

e(k 2 ) 

FIG. 3: Additional contributions to the Ai (qqi — > 77) amplitude. 




asM = 2m g . The resulting amplitudes are shown next: 



e(*i) 




Ai{xt, x 6 ,x u ) = ^7r4V3M ae 2|^ 9l (0)|e Al (fc 1 )e,(fc 2 )e p (P){(fc 2 -A :i )'' (1 - 4x 4 ) - 2^^^ 



+ (1 - 4a;4 - 2a; 5 ) 
+ (4sn 



2<? MP (M)" - 2^ 



M 2 



(hey [sT" + lkl J 2kl ) + (k 2 ey (> + M?_M-> 

which turns out to be independent of xq, xs,xio, due to the external momenta being on-shcll, and 
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A 2 {a) 



-ihir^£^a\*q qi (0)| e M (fci)e,(fc 2 )e p (P) 



e q sin 20\yK 1 ~ fl + 



M 3 / 2 

- J /i^8V3Ma|*g 9l (0)|e AJ (fc 1 )e,(fc 2 )e p (P) 

+2 9 ^ ( {kl er - 2k^-^l) + 2 9 ^ ( (k 2 er + 214 {kl0k2) 



M 



c q v K z . 



ena;P,-fa,-fe] 



M 



M 2 



M 2 



LCI,!' l.Phv\ / h p hP IcPh^ 



M 2 



+ (Mr 5 mp + 



M 2 



(a-1) 

}• 



(32) 



The coupling constants appearing in the above amplitudes are evaluated at the Mz scale [4, 20]. Here P = k\ + k 2 , 
(kid) 11 = ki V Q v>l and k\9k 2 = k\ ^O^ 1 " k 2 v , while M and P are the mass and the total momentum of the discussed 
quarkonium state, respectively. 

Both of the above amplitudes satisfy separately the usual Ward identities that follow from gauge invariance and the 
pole structure of the diagrams. Note that this would have been a problematic issue had we included in our deformed 
fermion lagrangian terms that contributed to the tree level fermion propagator. It is also worth noticing that, despite 
the dependence of the Feynman rules of eqs. (24) and (25) on the parameters xq,x$,xio -where x$ and x\o were 
associated to C, T violating contributions in the strong and hypercharge sectors -this dependence disappears in the 
amplitude Ai after evaluating the external momenta on-shell. In fact, the amputated Green function evaluated at 
generic momenta can be seen to depend on xq, x%, x\q, thus not contradicting the expectations of ref. [19], which 
stated the necessity of adding a term involving ^ vp to get one-loop renormalisability. Indeed, despite not contributing 
to the calculated amplitude, the associated terms in the basis of eq. (12) might influence other S-matrix elements or 
may be relevant to renormalisability properties. It is also worth to notice the dependence of the amplitude A\ in 
eq. (31) on the parameter in, which induces CPT violations in the weak sector. However, as can be seen from the 
diagrams of figure 2, the noncommutative vertices with fermion fields appearing in them only involve photons — no 
weak fields — , so that our computation is not affect by CPT violating interactions. 

In ref. [12], when computing the amplitudes for the above diagrams in the nmNCSM without the a-deformation 
of eq. (1) and without the additional deformations, Xj, of the fermionic terms, (22), it was found that the on-shell, 
amputated A\ amplitude, corresponding to the sum of the five quark-exchange diagrams in Fig. 2, was proportional 
to the on-shell, amputated amplitude A 2 of the boson-exchange Green functions represented in fig. 3, which could 
hint at a possible symmetry. This was made manifest by the fact that the S-matrix amplitudes in ref. [12] are given 
by some factor squared. In the present calculation, for a = 3, which is the one making the gauge sector of deformed 
nmNCSM one-loop renormalisable, this result cannot be recovered for any value of the free deformation parameters 
X4, X5, xq, x$, xio, in; in fact proportionality is only achieved for a = 1 and Xi = 0. The fact that we eliminated some 
of the freedom parameters X{ should not be relevant to this result, since we are concerned with on-shell, amputated 
amplitudes which are not affected by the field redefinitions that allowed us to eliminate some of the Xi. 

With respect to the possibility of considering terms deforming the tree- level fermion propagators, which have been 
ignored here, we would like to point out again that they would not only modify the on-shell conditions, but also 
cast serious doubts about the validity of the usual procedure of constructing the S-matrix amplitudes in terms of 
amputated diagrams; further considerations on these issues are out of the scope of this paper. 



Quarkonia decay rate 



To obtain the quarkonia decay rates we have to compute 

d 3 h d 3 k 

.. 'TZ2 I 



where M. is the total amplitude, i.e. A4 = A\ + A 2 . Summing over the initial spins and averaging over their final 
values in the square absolute value of the amplitude gives 
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3e 4 M 3 |5f 



(M 2 -M 2 ) 2 



A(e aP 9^) + B(8 ap 0P pkfkg) + C(8 ap 9? & k^k{ + ap 0" ^fcf ) + D(9 a/3 k?k$) 2 ] ; (34) 



i4 = 
S = 



M 2 



2C v K llZ M 2 + e 2 (4x 4 + 4^5 - 1) (M 2 - M|) + 2e 9 sin (26 w ) K~ nl (M 2 - M|) 2 , (35) 

(5a 2 - 22a + 25) C 2 K 2 7Z M 4 

+ 2e 2 q C v K^ z (M 2 - M§) M 2 ( - 28x4 - 28x 5 + a(12x 4 + 12x 5 - 4xu - 3) + 12xn + 7) 
+ e 2 sin 2 (20 w ) K 2 ^ (M 2 - M 2 ) 2 (5a 2 - 22a + 25) 

+ 2e 4 (l6x 2 + 8(4x 5 - 2x n - l)x A + 18a; 2 + 2x\ 1 + 4x n - 4x 5 (3xn + 2) + l) (M 2 - M§) 2 
- 2 sin (26V) if 777 (M 2 - Mf) (e 9 (5a 2 - 22a + 25) C v K llZ M 2 

+ - 28x4 - 28x 5 + a(12x 4 + 12x 5 - 4xn - 3) + 12x n + 7^ (Af 2 - Ml) ) , (36) 



C 



M 2 



(3a 2 -10a+ 11) C 2 K 2 jZ M 4 

+ 2e 2 q C v K llZ (M 2 - M 2 Z ) M 2 ( - 12x 4 - llx 5 + a(4x 4 + 5x 5 - 3x n - 1) + 5x n + 
+ e 2 sin 2 (26V) K 2 J7 (M 2 - M 2 Z ) 2 (3a 2 - 10a + ll) 

+ e 4 (M 2 - M 2 Z ) 2 (l6x 2 + 8(3x 5 - x u - 1)^4 + 12a; 2 + 4x 2 x + 2x u - 2x 5 (4xn + 3) + l) 

- 2 sin (26V) ^777 (M 2 - Mf) (e q (3a 2 - 10a + ll) C v K llZ M 2 

+ e g ( - 12x i -IIX5 + «(4x 4 + 5x 5 - 3xn - 1) + 5xu + 3j (M 2 - M%) ) , 



D 
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(3a 2 - 14a + 15) C 2 V K 2 ^ Z M A 

+ Ae 2 q C v K llZ (M 2 - Ml) M 2 ( - 82:4 - 8x 5 + a(4x 4 + 4x 5 - 2x u - 1) + 6x11 + 2) 
+ e 2 sin 2 (26V) ^ 777 (M 2 - M 2 ) 2 (3a 2 - 14a + 15) 



et(M 2 



Ml 



) 2 (16x1 + 8(4x 5 - 4xii - l)x 4 + 20x1 + 4x 2 n + 8x n - 8x 5 (3xn + 1) + l) 



2 sin (26 w ) if 777 (M 2 - Mf) (e 9 (3a 2 - 14a + Vb)C v K llZ 



M 



+ 2e q ( - 8x4 - 8x5 + a(4x 4 + 4x 5 - 2xn - 1) + 6xn + 2) (M 2 
Considering the quarkonia at rest and using the phase space integrals 



Mf)) 



M 2 ) 



d 3 hd 3 k 2 x4 . 



2 El 2E 2 - S ^ P - h -^ = A^-^2 



/ 



M 2 



^''Vfp-fa-fcXMfe) 2 = " - 



NC 

M 4 



2i^i 2i? 2 ii NC 
starting from (33), for general a, Xj, we obtain the following decay rate 

7T a 2 M 2 |*| 2 



\' ( ^ } 24' 



(37) 



(38) 



(39) 



nnmNCSM 



(gqi -> 77) = 



24 (M 2 - Mf ) 2 A 4 



NC 



(24A + 4M 2 (B + C) + M 4 D S jE 2 - (24A + 2M 2 (B - 2C)j B 2 .(40) 
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The above rate turns into eq.(12) from [12] for a = 1, xi = 0. We would like to analyse eq. (40) and study the effects 
of having added to the action of the nmNCSM the extra terms that have been discussed in the previous sections, 
in order to compare the results with those calculated in rcf. [12]. In order to obtain some numerical values, we will 
look for the maxima and minima of the decay rates in the allowed region for K Zll and i"T 777 -see figure 1- with the 
assumption that the Xi parameters are "natural" and only take values between zero and one, and with two possible 
scenarios for the dimensionless constants Eg and Bg : either both of them arc of order one (space-time and space-space 
noncommutativity), or Eg = and Bg is of order one (only space-space noncommutativity). Moreover, as was done 
in ref. [12], we will consider that the scale of noncommutativity varies between Anc = 0.25 TeV and Anc = 1 TeV. 

In order to compute branching ratios, we use the following data taken from [20]: in the J/ip case, r cxp (J/^ — > 
e+e") = (5.55 ± 0.14 ± 0.02) keV and T^ t p (T) = (93.2 ± 2.1) keV, whereas for the T case, r ox P (T -> e+e") = 
(1.340 ± 0.018) keV and r°* p '(T) = (54.02 ± 1.25) keV. Recall that the wavefunction at the origin (0) is related to 
the lepton decay rate by eq. (30). 

First, since in the minimal NCSM (mNCSM) there are no .Z77 an d 777 couplings, the formula (40) can be used to 
recover the branching rations in the minimal NCSM -see ref. [12], eqs. (19,20)- by taking K Zlll = K 111 = 0, which 
yields, for Ej = B 2 = 1 and 0.25 TeV < A NC < 1 TeV, 



5.1 • 10- 13 < BRfflCJ^ < 1.3 • 10" 
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4.6 • 10- 12 < BRftl° 7 T < 1.2 • 10- 



(41) 



The values at Eg = 0, Bg = 1 are suppressed by a factor of 3/10. 

In ref. [12], the computation in undeformed nmNCSM was also done for Eg = Bg — 1, yielding 



T a 1 (qq 1 -> 77) 
T(qqi -» £+1-) 



24 



M 



A 



NC 



■ sin 29w K-y 



\Mz) 



CyJKz'yj 



(42) 



This corresponds to setting a = 1, Xi = in the deformed version of the nmNCSM of this paper. 
Computing the maximal values of the above rate for 0.25 TeV < Anc < 1 TeV yields 



3.1 • 10~ 12 < BR, a u]' x<= 1 ° < 7.8 • 10- 10 , 

~ [J/y— »'yy\,max ~ ' 



1.7 • 10- 11 < BR " =1 - x r° < 4.3 • 10" 9 , 

~ [T— >77j .max ~ ' 



(43) 



as was obtained in eqs. (23,24) from [12]. 

On the other side, the minimal values in the same range of Anc are 



2.3 • 10~ 13 < BR, a =\ x >=° . <5.9-10- n , 1.0 • 10" 26 < BR ° =1 ' x ' =0 . < 2.6 • 10" 24 . 

~ [J/ip— >77j,mTO ~ ' ~ [T— >77j,mm ~ 



(44) 



By taking Eg — 0, Bg — 1 all the results are suppressed by a factor of 3/10. 

To see the effect of adding to the action the extra pure gauge term associated with the choice a = 3 in eq. (5), 
corresponding to the deformed nmNCSM with renormalisable pure gauge interactions, we can start from (40) and fix 
a = 3 and X4 — x 5 — xn — 0. The resulting expression is more complicated than (42), 



r a=3 (g qi 



77) 



M 4 



T(qqi -> l+l 

R = 3(M 2 - IV, y 



48e^ c (Mi-Mlr iRBd+SEe 
1 M y ) 2 e„ i + 36sm 2 (2e w )K^(M 2 
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2 e 2 + 20c q v M 2 K z ^ (M 2 



M 2 Z ) e 2 
4^2 



4 sin (29 w ) K„ nl (M 2 - M%) (5 (M 2 - M 2 Z ) e 2 + 18c q v M 2 K z ^) e q + 2,&c q v z M A K 2 Zll 



S = 7 (M 2 - M\ 



+ 36sin 2 (2<9 w )iC 2 [M 



M\) % 2 



20c q v M 2 K 



z-rr 



(M 2 



M z z ) 

4^2 



+ 4 sin (28 w ) i^ 777 (M 2 - M|) (he 2 (M 2 - M 2 Z ) - 18c q v M 2 K z ^) e q + 36c q v ' M 4 K 2 Z ^ 



(45) 



The maximal values in the allowed regions for K Zll and K 111 of the branching ratios, for 0.25 TeV < Anc < 1 TeV, 
El — B 2 — 1 become 



2.4 • 10~ 12 < BR, °= 3 ' x '=° < 6.3 • 10- 10 , 



7.5 • 10~ n < BR, a ~ z - x '-° < 1.9 • 10~ 8 . 

~ [T— >77j, max ~ 



(46) 



whereas for the minimal values we get 



5.1 • 10~ 13 < BR, a = 3 ' Xi =° . < 1.3 • 10- ln , 



4.6 • 10~ 12 < BR,^ 3 ' x r° < 1.2 • 10~ 9 . 

~ [T—>77j, min ~ 



(47) 
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They remain more or less in the same order of magnitude as those of eq. (43), safe for the minimal values of the T 
branching ratio, which are hugely increased. In particular, the minimal values are always above those of the mNCSM 
of eq. (41). In the case Eg = 0, B 2 = 1 the rates are again suppressed but in a different way: 

1.9 • 1(T 12 < BR, a = 3 ' x '=° < 4.8 • 1(T 10 , 2.3 • l(T n < BR [7 T 3 < T° < 5.9 • 1(T 9 , 
3.5 • 10~ 14 < BR, ^ 3, . <9.0-10~ 12 , 1.0 • 10" 13 < BR " =3, x ; = ° < 2.6 • 10~ n . (48) 

In particular, the minimal values arc no longer above the ones corresponding to the mNCSM with Eg = 0,Bg = 1, 
but there is still a dramatic increase in the minimal possible values for the T branching ratio. 

So far, the effect of the extra pure gauge term in the deformed nmNCSM with respect to the undeformed version 
is essentially an enhancement of the minimal allowed values of the T branching ratio. 

Next we move on to see the effect of the extra terms in the fermionic action depending on X4,X5,Xn, when these 
parameters take natural values between -1 and 1. We set to calculate the maximal and minimal values that the 
branching ratios obtained from eq. (40) can have for Xi € { — 1, 1} within the allowed values of K Zll and K 111 . 

In the case Eq = Bg = 1, we obtain the following results, as before for scales 0.25 TeV < Anc < 1 TeV: 

4.5 • 10~ n < BR^ 3 ' < 1.2 • 10~ 8 , 5.6 • lO" 10 < BR, a = 3 > < 1.4 • 10~ 7 , 

a=3, „ Bi2 a=3, _ g (4g) 

[J/ip— >u\,min 1 [T— >77j,7mn V ' 

As an explanation of the previous values, for example the first expression 4.5 • 10~ n < B R^jJ^'J^^max < 1.2 • 10~ 8 
means that the maximum of the J/ip branching ratio, for the range of scales 0.25 TeV < Anc < 1 TeV, for all 
Xi E { — 1, 1} and for all allowed K Zll and -^777, varies between 4.5 • 10~ n -reached at the scale Anc = 1 TeV, and 
1.2 • 10 -8 , reached at Anc = 0.25 TeV. Note that the dependency of the rates on Anc factorises, as seen in eq. (40). 
Each value is reached for a particular value of Xj, -Kz 77 , if 777 : 

Maxima (50) 
J/ip : x 4 = -1.00, x 5 = -1.00, xn = 1.00, K Zll = -0.254, AT 777 = 0.129, 
T : x 4 = -1.00, x 5 = -1.00, an = 1.00, K Zll = 0.00950, AT 777 = -0.576, 

whereas for the minima the numerical results are not reliable due to precision issues. Figures 4 through 7 show 
the resulting branching ratios as functions of Kz^-y and iif 777 for the particular values of the Xj that yielded the 
above maxima and (approximate) minima, respectively, at the scale Anc — 1 TeV. The maxima of figures 4 and 
5 correspond to the values appearing at the left of the inequalities in eq. (49). The minima of figures 6 and 7 are 
practically zero. 

For completeness we show the corresponding results in the Eg = 0, Bg — 1 case: 

2.0 • 10- 11 < BR, a = 3 ' < 5.1 • 10- 9 , 2.0 • 10- 10 < BR "= 3 ' < 5.0 • 10" 8 , 

~ [J/ip— *77j, maa; ~ ' ~ [T— »-77j ,max ~ ' 

BR a=3, |x ( |~i ^ BR^ 3 - 1 *^ 1 ~0, (51) 
The previous values are reached for the following values of the parameters: 

Maxima (52) 
J/ip : xi = -1.00, x 5 = -1.00, an = -1.00,^77 = 0.00950, X 777 = -0.576, 
T : x 4 = 1.00, x 5 = 1.00, xn = 1.00, K Zll = 0.00950, if 777 = -0.576, 

where again we do not display the position of the minima since our numerical results are not reliable. The effect of 
the Xi terms is clearly to allow for much larger maximum values and much lower minimum values of the branching 
ratios with respect to the results in the undeformed nmNCSM or the mNCSM. The maximum values are increased 
up to 2 orders of magnitude with respect to the mNCSM result of eq. (41) and one order of magnitude with respect 
to the results in eq. (43) corresponding to the undeformed version of the nmNCSM. 

Despite the fact that the allowed values for the minima experience an important decrease, for typical values of 
the parameters K Zll , K 711 , x 4 , x 5 , xn the rates are enhanced with respect to the mNCSM result. To justify this 
claim we have computed the branching ratios for random values of the above parameters in the allowed region for 
K Zl ~ 1 , K 171 displayed in figure 1 and for Xj € { — 1,1}. For up to one million configurations of these parameters, at 
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FIG. 4: J/'i! —* 77 branching ratio as a function of -Kz 77 and -K" 777 , for Eg,Bg ~ 1, #4 = X5 = — asu = — 1, at the scale 
A NC = 1 TeV 




FIG. 5: T — » 77 branching ratio as a function of JCz 77 and if 777 , for Eg,Bg ~ 1, 0:4 = 2:5 = — £11 = —1, at the scale 
A NC = 1 TeV 



the scale A NC = 1 TeV, we found that, for Eg = Bg = 1, (E$ = 0, Bj = 1), 96(88)% of the configurations yield J/tp 
branching ratios larger than their corresponding values in the mNCSM, whereas in the T case the percentages are 
97(89)%. The percentage of configurations yielding a 10 x increase over the mNCSM values are 44(42)% in the J/tp 
case and 55(40)% in the T case. There is a 50x increase for 2(4)% of the J/tp configurations and for 4(3)% of the T 
configurations. 

We have also estimated the portion of the parameter space which, at the scale Anc = 1 TeV, yields values of the 
branching ratios which are larger than the maximum values in the undeformed version of the nmNCSM, which are 
given in eq. (43). The percentage of configurations that satisfy this requirement for Eg = Bg = 1, (Eg =0, Bg= 1), 
is given in the J/tp case by 62(55)%, whereas in the T case it is equal to 84(67)%. 

To complete our numerical estimates, we should consider the case when the parameters Xi take values that are 
necessarily non-natural, i.e, not necessarily restricted to be between -1 and 1. An analysis of the dependence of the 
branching ratios on the parameters i^z 77 , if 777 , xi shows that, independently of the Xi, the branching ratios as a 
function of iCz 77 and -ftf 777 define a concave parabolic surface (in the J /\& case) or a hyperbolic surface (T case). 
Thus the maxima always appear in the boundary of the region of the parameter space, as has happened with our 
results in eqs. (50) and (52). It is clear that allowing a wider range of the Xi will directly yield greater maxima; since 
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FIG. 6: J 1^1 — > 77 branching ratio as a function of iCz 77 and K lnl , for Eg, By ~ 1, an = 0.25, = —2.85 ■ 10 7 ,m = 
2.86 • 10~ 7 , at the scale Anc = 1 TeV. The horizontal plane represents the mNCSM branching ratio of 5.1 • 10 -13 . 




FIG. 7: T — > 77 branching ratio as a function of Kzj-y and K~ nj , for Eg, Bg ~ 1, X4 = 0.25, x$ — —2.85 -10 7 , xn = 2.86 ■ 10 7 , 
at the scale Anc = 1 TeV. The horizontal plane represents the mNCSM branching ratio of 4.6 ■ 10~ 12 . 

the dependence of the branching ratios is quadratic in the Xj, we expect that an increase of an order of magnitude 
in the range of the Xi yields a two orders of magnitude increase in the maximum values of the branching ratios. By 
doing some numerical computations this seems to be roughly the case (modulo a factor between 0.45 and 0.77). In 
particular, taking xi € {-100, 100}, we get, for Eg = Bg = 1 and for scales 0.25 TeV < A NC < 1 TeV: 

3.5 • 10~ 7 < BR, ,7 3 ' Nl r 100 < 9.4 • 10^ 5 , 3.2 • 10~ 6 < BR [7 T 3 < l * i| ~ 100 < 8.1 • 10~ 4 . (53) 

~ [J ftp— ry"f\,max ~ 5 ~ l±— >77j ,max ~ V" i 

The values at the scale A = 0.25 TeV are within the experimental bounds: BR(J/ip — ► 77) < 2.2 x 10~ 5 [20] and 
Bi?(T^ 77 ) < 10- 4 [21]. 

DISCUSSION AND CONCLUSIONS 

In this paper we have introduced a deformed fermion lagrangian for the nmNCSM and we have applied it to the 
phenomenological estimate of the quarkonia decay rates into two photons. 
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First of all, let us recall that the non-minimal version of the NCSM gauge sector includes traces over the represen- 
tations of all the massive particle multiplcts with different quantum numbers that appear in the total lagrangian of 

the model with covariant derivatives acting on them, i.e., with terms of the type i'tp * fl * %/j , (D^Qy ★ (1)^$), [4, 5]. 

Second, we have argued why the use of the nmNCSM, and in particular its deformed version, should be favoured 
over the mNCSM; this is because of the proven renormalisability of the pure gauge sector and the existing hints 
that this property might not be spoilt when the effects of matter loops are taken into account. The motivations 
for introducing deformation terms in the fermion lagrangian come essentially from our lack of knowledge of the 
renormalisability properties of the matter sector of the full nmNCSM, and also from the fact that, for the gauge 
interactions, renormalisability -and even finiteness for the first noncommutative corrections- was only achieved by 
adding a deformation term to the starting lagrangian. Thus, our deformed fermion lagrangian should be taken as 
an effective lagrangian for the nmNCSM as long as the renormalisability properties are unknown. From the effective 
theory point of view, the values of the Xi parameters could be constrained by experimental measurements. 

It should be noted that the photon polarisation is known to be modified by noncommutativity, which causes vacuum 
birefringence [22], at least in the standard approach to noncommutative theories. This follows from computations 
of the one-loop photon self-energy. It should be interesting to analyse this issue in the enveloping algebra approach. 
Of course our new fermion-photon interactions, coming from the x^-dependent terms, will also affect the photon 
polarisation at one-loop, which could be used for further experimental tests. However this is not straightforward 
since, according to rcf. [23], the photon two point function, due to gauge invariance, will only be modified at order 9 2 
and beyond, so that a consistent computation would imply the use of SW maps up to 0{6 2 ). 

Aside from possible experimental measurements, there is hope that further investigations about renormalisability 
properties could impose constraints on the free deformation parameters Xi; it could happen that some or all of them 
were fixed uniquely, as happened for the gauge sector of the nmNCSM, and in that case the new fermion sector would 
not be effective but part of a theory well defined in the ultraviolet. 

In deriving the deformed fermion lagrangian, we did not consider terms that altered the tree-level fermion propaga- 
tor. This was done because, on the one hand, we should expect noncommutative effects to appear as weak quantum 
corrections, and thus we do not find desirable to break the usual Lorentz-invariance matter dispersion relation p 2 = m 2 
at tree level. On the other hand, we wanted to apply the framework to the computation of decay rates, which implies 
the calculation of S-matrix elements. In ordinary space-time one uses the LSZ formula, which relies on general prop- 
erties of the pole structure of the Green functions of the theory (and in particular the 2-point function) which follow 
only from Poincare invariance. In noncommutative space-time the usual Lorentz invariance does not hold, and thus a 
proper all-order definition of S-matrix elements seems challenging (in the case of noncommutative theories that do not 
make use of Seiberg-Witten maps, some advances have been done in this respect, see for example ref. [18]). Neverthe- 
less, since we were computing S-matrix elements at tree level and we were not altering the tree-level 2-point function, 
we believe that the usual LSZ formalism should hold at this level. As a consistency check, our S-matrix amplitudes 
satisfy the usual Ward identities associated with U(l) gauge invariance, which are derived in ordinary space-time as a 
consequence of gauge symmetry and the pole structure of the diagrams. Nevertheless, a deeper understanding of the 
S-matrix and the LSZ formalism in noncommutative theories is still needed. We also recall that some of the terms 
considered in eq. (12) violate CPT exclusively in the weak sector; also, there appear C,T violations in the strong and 
hypercharge sectors. This could be of phenomenological interest for searches of Physics beyond the Standard Model; 
note that the violations are very small since they appear purely as noncommutative effects. 

Concerning our results for the quarkonia decay amplitudes, it should be noticed first that the on-shell amplitudes 
turned out to be independent of the C, T violating terms t s ,ti ; also, though they are dependent on in, which is 
associated to CPT violations in the weak sector, no CPT violating interactions contributed to the result. Ref. [19] 
argued that to get one-loop renormalisability in noncommutative QED, the term of cq. (12) should be added to 
the bare lagrangian. This does not conflict with the fact that the quarkonia decay amplitudes are not apparently 
influenced by the terms with three 7^ matrices in eq. (12): first, we absorbed the contributions of ty in the terms 
£4 — t$, ts, t\o, tu by using field redefinitions, and, though the on-shell amplitudes did not depend on the parameters 
X6, x 8 , X10, the amplitudes evaluated at arbitrary momenta did. With respet to the numerical results obtained for the 
decay rates, we showed that the effect of the extra gauge term in the action of the extended nmNCSM is essentially 
to raise the minimum allowed values of the T — > 77 decay rate, masking the destructive contributions that were 
found in the undeformed nmNCSM in rcf. [12]. However, when the new terms of our deformed fermion lagrangian are 
taken into account, this effect disappears, and both the minimum and maximum allowed values for the decay rates 
experience high decreases and increases, respectively. 

Having very low allowed values for the quarkonia decay rates is no good news for the possibility of comparison 
with experiments in order to confirm or falsify the theoretical predictions: in principle, non-zero qqi — ► 77 decay 
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rates could be taken as a signal of noncommutativity, but if the models allow for extremely small values it would be 
difficult to discard them. However, the panorama is more promising because, despite the minimum possible values 
are very small, for most configurations of the parameters (between 88 and 97 % at A NC = 1 TeV for natural values 
of the parameters) we obtain branching ratios that are greater than the ones that were computed for the mNCSM. 
Furthermore, for a big portion of the parameter space (between 55 and 84 % under the same conditions), the decay 
rates, computed from deformed nmNCSM, are actually larger than the maximum values that were found in the case 
of the undeformed nmNCSM. 

Thus, in general we get constructive contributions to the decay rates and the model allows for greater values of 
the branching ratios that the ones that had been previously found; the maximum values increase by up to two orders 
of magnitude. Possible future studies of the renormalisability of the matter sector could help to restrict the allowed 
values in the parameter space, as happens in the gauge sector with the ambiguity parameter a of eq. (5) forced to be 
equal to 3, and this could make the model more predictive or falsifiable. Again, we recall the result of ref. [19], which 
in the QED case argues that one-loop renormalisability demands to add only the term tr of eq. (12) to the lagrangian; 
this may also happen in the NCSM but the result cannot be directly extrapolated since, after the expansion with 
the SW map, the NCSM lagrangian is not given by a sum of lagrangians for the different gauge groups due to the 
appearance of interactions between the different gauge fields. 

Todays existing experimental limit for the branching ratio of the J/ip — > 77 decay can be found in "Review of 
particle physics," under the C symmetry violating modes, and is BR(J/ip — > 77) < 2.2 x 1CP 5 [20]. 

With respect to the T case, as it was commented in rcf. [12], the existing limit for the branching ratio of the 
T — ► 77 decay comes from a very old CLEO-III experiment [21] and it indicates that with present data the detection 
of BR(T — > 77) below 10~ 4 would be hopeless. 

As follows from eq. (53), the previous experimental bounds are reachable in our model at the scale A = 0.25TeV 
for values of Xi tuned to be around 100, and for lower scales for higher values of the X{. Thus there is some hope that 
the phenomenology of our model could be relevant. 

Furthermore, despite the fact that all experiments including hadrons are extremely hard to perform and analyse 
due to the huge background signals, the large number of heavy quark-antiquark pairs harvested at LHCb; i.e. 10 12 
BB pairs per year [24], and probably (10 14 DD; 10 18 KK), give us hope that experimental branching ratios 
BR(T — > 77) <~ 10 -9 and BR(J '/ip — > 77) ~ 10~ n could be accessible, thus reaching our maximum predicted values 
for the rates for \xi\ <~ 1 (46-51). Certainly these experiments would produce at least reliable (and much lower) 
bounds; perhaps they could even measure these processes, depending on the scale of noncommutativity. We hope 
that the importance of a possible discovery of space-time non-commutativity will convince experimentalists to look 
for SM forbidden decays in hadronic physics. 
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